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We propose a new approach to the Kondo lattice in order to describe simultaneously the Kondo 
effect and the local magnetism. This approach relies on an original representation of the S — 1/2 
impurity spin in which the different degrees of freedom are represented by fermionic as well as bosonic 
variables. The magnetic instability is found at Jc/D = 0.67 in definite improvement compared to 
usual fermionic mean-field theories. The dynamical susceptibility contains an inelastic peak in 
addition to the standard Lindhard-type Fermi liquid contribution. The Fermi surface sum rule 
includes nc + 1 states corresponding to large Fermi surfaces. 



An outstanding feature of heavy-Fermion systems is 
the coexistence of Fermi liquid-type excitations with lo- 
cal magnetism resulting from RKKY interactions among 
spins. For instance, the dynamic magnetic structure fac- 
tor measured in CeCug and CeRu2Si2 by inelastic neu- 
tron scattering (INS) indicates two contributions [0J|]: a 
q-independent single-site slow component typical of the 
localized excitations of Kondo-type, superimposed on a 
strongly q-dependent intersite fast component, reflect- 
ing the magnetic correlations due to RKKY interactions. 
The former is peaked around uj = (quasi-elastic peak) 
while the latter is peaked around ojq 7^ (inelastic peak) . 
Another important probe is provided by the de Haas- 
van Alphen (dHvA) experiments. The results ||] ob- 
tained by different groups on various compounds indicate 
heavy effective mass. They also agree as to the existence 
of large Fermi surfaces in the magnetically-disordered 
phase. Even though the charge degrees of freedom are 
frozen, the localized electrons seem to contribute to the 
Fermi surface sum rule together with the conduction elec- 
trons. 

This coexistence of Fermi liquid-type excitations with 
low energy magnetic fluctuations is likely to stem from 
the nature of the screening of the localized moments in 
the Kondo lattice. The interesting domain is reached 
when I <Tk <T3 {I is the order of the RKKY interac- 
tion, Tpc the Kondo temperature and T3 the temperature- 
scale at which the number of thermal states for con- 
duction electrons just equals the number of impurities). 
Then, at Tk the available conduction electrons are ex- 
hausted Q before achieving complete screening (incom- 
plete Kondo effect) leaving residual unscreened spin de- 
grees of freedom on the impurities. 

Traditionally, the spin is described either in fermionic 
or bosonic representation. If the former representation, 
used for instance in the 1/iV expansion of the periodic 
Anderson-PAM [|| or the Kondo lattice-KLM |§] mod- 
els, appears to be well adapted for the description of the 
Kondo effect, it is also clear that the bosonic representa- 



tion lends itself better to the study of local magnetism. 
Quite obviously the physics of heavy-Fermions is domi- 
nated by the duality between Kondo effect and localized 
moments. This constitutes our motivation to introduce a 
new approach to the KLM which relies on an original rep- 
resentation of the impurity spin 1 /2 in which the different 
degrees of freedom are represented by fermionic as well 
as bosonic variables. The former are believed to describe 
the Fermi liquid excitations while the latter account for 
the residual spin degrees of freedom. 

Let us list the main results obtained in this paper: 
(i) the Fermi surface sum rule includes ric + 1 states 
which means that the Fermi surface volume includes a 
contribution of one state per localized spin in agreement 
with dHvA experiments, (ii) the dynamical susceptibil- 
ity contains an inelastic peak in addition to the standard 
Lindhard-type Fermi liquid contribution as observed in 
INS, (iii) the instability towards long-range magnetic or- 
der appears at Jc/D — 0.67 significantly reduced to the 
value Jc/D = 1.0 obtained in purely fermionic mean- field 
theories. While (i) has been already obtained in previous 
work, we stress on the fact that (ii) and (iii) constitutes 
new results established by this supersymmetric approach. 

We propose the following mixed fermionic-bosonic rep- 
resentation of the spin S = 1/2. The two states constitut- 
ing the basis may be written: 1 1/2, 1/2} (Xfe|+r/|)|0} 

and |l/2, -1/2) = (Afej + F/|)|0), where and /t are 
respectively bosonic and fermionic creation operators, X 
and Y are parameters controlling the weight of both rep- 
resentations and |0) represents the vacuum of particles: 
bcr |0) = fa |0) = 0. The spin operators are then given by 

= E..' blr-^'ba^ + nT-^,U' = + Sf, where 
{a ~ (+, — , z)) are Pauli matrices. Sf, and S/ are equiv- 
alent to Schwinger boson and Abrikosov pseudo-fermion 
representations of the spin respectively. 

One can easily check that the representation that we 
propose satisfies the standard rules of SU{2) algebra: 
1/2, ±1/2) are eigenvectors of 5^ and with eigen- 
values 3/4 and ±1/2, [S+,S-] = 2S'' and [S'-,S^] ^ 
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±5* provided that the following local constraint is sat- 
isfied: X](T ^CT^'T + /it/o" ~ nb + Uf = 1, which implies 
+Y'^ — 1. In order to eliminate the unphysical states 
as (^Xb'l + |0), we need to introduce a second lo- 

cal constraint and take: Q = J^a ^{fl^-iy + ^-afcr) = 
Let us consider the three-dimensional KLM near half- 
filling [ric < 1). The hamiltonian is: 



(1) 



In the representation introduced before, the partition 
function can be written as the following path integral: 



(2) 



with Co = E»a(cl<T5rQ<x + fladrha + h\^drh„) 

and W = Efe<T efc4<TCfc<T 



where J(> 0) is the Kondo interaction, and two time- 
independent Lagrange multipliers e/i and Ki are intro- 
duced to enforce the local constraints. 

Performing a Hubbard- Stratonovich transformation 
and neglecting the space and time dependence of the 
fields in a self-consistent saddle-point approximation, we 
have : 

Z = / d-qdT^*C (ctq, Aq, 77, 77*) Z(77, 77*) 
Z(77, ??*)=/ Vc,„Vf,„Vb,„ exp [- dT {Co + H') 



Hq being hermitian, the matrix transforming the 
original basis ijj''' = (pc'b^) to the basis of eigenvectors 
$t = (at/jt^t) is unitary {UW = C/tj/ = 1). $t ^ 
^jj^U^ with a supersymmetric matrix, and are 
the fermionic eigenvectors whose eigenvalues, determined 
from det [{a - E) - a{b - E)-'^ p] = 0, are: 



E^^ 



(£fc + £/) T V(£fc-£/)'+4K+W*) 



7^ is the bosonic eigenvector whose eigenvalue, deter- 
mined § from det [{b ~ E) - p{a - Ey'^a] = is = 

In the scheme we propose, ctq and Aq are slow vari- 
ables that we determine by solving saddle-point equa- 
tions, while 77, 77* are fast variables defined by a local 
approximation. As we will see, the latter approximation 
incorporates part of the fluctuation effects. Indeed, per- 
forming the functional integration of (|^) over the fermion 
and boson fields [Q yields a superdeterminant (SDet) 
form written as follows: 



Z(77 , 77*) = SDet{dr + H) , 



where SDet{dr + H) 



Det{G-^ - aDp) 
Det{D~^) 



(4) 



=dr + a and D'^ =dr + b 



Expanding to second order in 7;, 77* allows us to define 
the propagator G^^* (k, iujn) associated to the Grassmann 
variable 77 and hence the closure relation for = {rjrj*): 



Xq 



4 G^')* (k, iw„) 



(5) 



with C (cto, Ao,77, 7;*) = exp[/3 ^Xq - Sjl^IUL - ^Y^^ri^ 

V ho 

(To ~ 
Hq^ \ (To Ek ri 
77* Sf 



J with Gijij. (k, ia;„) = 



J 



(3) 



Note the presence of a Grassmannian coupling rj between 
Cicr and bia- , in addition to the usual coupling ctq between 
Cia and fia responsible for the Kondo effect. Ho is of the 
a a 



type ^ ^ J in which a, 6 (p,(t) are matrices consisting 

of commuting (anticommuting) variables. Note the su- 
persymmetric structure of the matrix Hq similar to the 
supermatrices appearing in the theory of disordered met- 
als 01 . The second constraint is automatically satisfied 
in this scheme since there is no off-diagonal term mixing 
t end J, states. 



and n2, = i^ 



[1- jn2,(k,i.;„)] 

Gcc(k + q, iujn + iuj^)D{q, 7lj„) . 



Contrary to ||8| which in the case of the underscreened 
Kondo impurity model, assumes Xq = leading to a two- 
fluid description, the closure equation (^) that we intro- 
duce defines a finite Xq. This parameter Xq plays a ma- 
jor role in controlling the relative weights of fermion and 
boson statistics. It is directly connected to the X and Y 
parameters introduced in the initial representation of the 



states : X"^ = ^^/{al + x^) and = (jl/{al + xl). 

The resolution of the saddle-point equations, keeping 
the number of particles conserved, leads to: 



VF = -i? exp [-1/(2 Jpo)] , 
2pq{(jI + xl) 



1 = 



-yp 

D 



(6) 
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where yp = fi — Sf and po = i/'2D is the bare density of 
states of conduction electrons. From this set of equations, 
we find : Sf — 0. 

The resulting spectrum of energies is schematized 
in Figure At zero temperature, only the low- 

est band a is filled with an enhancement of the den- 
sity of states at the Fermi level (and hence of the 
mass) unchanged from the standard expansions: 



piEp) 



= 1 



K + 4) 



= 1 



D 



> 1. This large 



Pa {-yp) 

mass enhancement is related to the flat part of the a 
band associated with the formation of the Abrikosov-Suhl 
resonance pinned at the Fermi level. While this feature 
was already present in the purely fermionic description, 
it is to be noted that the formation of a dispersionless 
bosonic band within the hybridization gap is an entirely 
new result of the theory. 
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FIG. 2. J/_D-dependence of the coupling — {"qri"} fixing 
the relative weight of fermion and boson statistics. The unit 
on the vertical-axis is . 




FIG. 1. Sketch of energy versus wave number k for the 
three bands a, /3, 7 resulting of the diagonalization of super- 
symmetric Ho- Note the presence of a bosonic band 7 sepa- 
rating the fermionic bands a and p. Note also the flatness of 
the band at fe = /cf reflecting large effective mass. 

The relative weight of boson and fermion statistics in 
the spin representation is related to Xq : ni,/nf = Xq/ctq. 
It is then interesting to follow the J-dependence of Xg. 
The result is reported in Fig. |2| This bell-shaped curve 
can be interpreted in the light of the exhaustion princi- 
ple mentioned in the introduction. In the limit of large 
J, the Kondo temperature-scale Tk = Dexp[—l/{2JpQ)] 
is of order of the bandwidth. One then expects a com- 
plete Kondo screening as can be checked by remarking 
that the weight of c in the a quasiparticle at the Fermi 
level (noted v^^ ) just equals the added weights of f and 
b at the Fermi level (respectively noted uf. and pf): 
+ Pi) = vll[ol + xg) = 1. The^Kondo ef- 
feet being complete in that limit, there is no residual 
unscreened spin degrees of freedom: it is then natural to 
derive a zero value of x\ (and hence of Uh). The opposite 
limit at small J corresponds to the free case of uncoupled 
impurity spins and conduction electrons. It also leads to: 
Xq = 0. The finite value of x\ between these two limits 
with a maximum reflects the incomplete Kondo screening 
effect in the Kondo lattice, the unscreened spin degrees 
of freedom being described by bosons. 



Largely discussed in the litterature ||Tl|] is the ques- 
tion concerning the Fermi surface sum rule: do the local- 
ized spins of the Kondo lattice contribute to the counting 
of states within the Fermi surface or do they not? De- 
pending on the answer, one expects large or small Fermi 
surfaces. The supersymmetric theory leads to a firm con- 
clusion in favour of the former. One can check that the 
number of states within the Fermi surface is just equal 
to ric -1- rib -|- , i.e. ric + 1. The Fermi surface volume 
includes a contribution of one atate per localized spin in 
addition to that of conduction electrons This con- 
clusion which appears reasonable if one recalls that the 
KLM is an effective hamiltonian derived from the PAM, 
has been reached before by several other authors ||ll|,|l2| . 
We do think that it is a good sign to recover previously 
established results if they are correct. 

The dynamic spin susceptibility x(q, iS) is calculated as 
usual from the second derivative of the generating func- 
tion lnZ{Hi) with respect to a magnetic source Hi. The 
corresponding spin-spin correlation function is given by: 
< {Sf{t) + Sb{t)){SfiO) + SbiO)) >. Projecting it in the 
eigenbasis {a, (3, 7} and omitting negligible terms, one 
gets two contributions to (q, cu) : 



x" (q, = x" aa (q, + x" a7 (q, w) 



(7) 



The first one is of Lindhard type and corresponds to 
electron-hole pair excitations in the a band: 



X"aa(q,t^) ~ piEp)—^ 

qv*p 



(8) 



where Vp is the renormalized Fermi velocity. The second 
one corresponding to the excitations from the a to the 7 
band has a peak at a frequency of order Tk- 



X"a^(q,w) Pq 



1 



w2(cro +Xo ) 



(^ - yp) (9) 



We believe that these two contributions may account for 
the quasielastic and inelatic peaks observed in INS. 
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FIG. 3. J/-D-dependence of the renormalized populations 
iVo and (where iVo = Nq.Tk and Wb = SN^.Tk) of the 
bosons in the condensate and the excited levels respectively. 
The intersection of the two curves determine the critical value 
Jc/D = 0.67 below which long-range ferromagnetic order ap- 
pears. Note the first order transition at Jc/D with A''o (and 
5Nb) jumping to (and from) zero at this value. 

One of the crucial questions raised by the Kondo lat- 
tice model is the possibility of long-range magnetic order 
at low J. So far, all the results have been obtained as- 
suming that no Bose condensation occurs. Let us now 
reconsider the problem allowing the bosons b condense. 
At the level of approximation that we use, bosons be- 
have as free particles; so at T 0, the bosons are ei- 
ther all in the excited levels or all in the condensate. In 
the latter case, we evaluate the number of bosons A^o 
in the condensate to No — exp[— l/2po>/]- This value 
should be put together with the number of bosons 5N\j 
in the excited levels in the absence of Bose condensa- 
tion: 5Nb = 2poxl/{-yF) = xl/{D^e^p[~l/2poJ]). It 
turns out that comparing the free energy of both states 
is equivalent to compare Nq and SNi,. Figure || repre- 
sents the J-dependence of the two populations. The in- 
tersection of the two curves defines the value of J below 
which the bosons condense. The ferromagnetic insta- 
bility associated to this Bose condensation appears at 
Jc/D = 0.67. This value is significantly reduced com- 
pared to the fermionic theory mean-field value of order 
1.0 The reduction of Jc/D that we predicted is 

similar to what is found in the related problem of the 
ID Kondo necklace jl^ for which numerical calculations 
also lead to smaller critical value {Jc/D = 0.382 and 
0.24 within real-space renormalization group and finite 
size scaling analysis respectively) compared to mean-field 
value {Jc/D = 1.0). It is obvious that the sharp tran- 
sition of iVo and SNi, at Jc associated with a first-order 
transition is an artifact of our approximations. Further 
refinements including fluctuations would introduce inter- 
actions among bosons, as in the Bogoliubov problem. We 
expect that it will produce progressive depiction of the 
condensate towards the excited levels transforming the 



ferromagnetic transition from first to second order. With 
a few exceptions, most of the heavy Fermions are not 
ferromagnetic but rather antiferromagnetic. The present 
calculation cannot directly apply. However, we believe 
that the discussion here is very useful since (i) it draws 
the general lines of the study which could as well be ex- 
tended to the case of antifcrromagnetism by considering 
a staggered instead of an uniform condensate and (ii) it 
can be considered as a test of the theory by comparing 
the prediction of the ferromagnetic instability with the 
results of numerical calculations. 

In summary, we have shown that the supersymmet- 
ric method proves to be a powerful tool to account for 
both Fermi-liquid excitations and residual spin degrees of 
freedom in heavy-Fermion systems. The key results are 
the prediction of a magnetic instability at Jc/D — 0.67 
significantly reduced compared to purely fermionic theo- 
ries, the presence of two contributions in the dynamical 
spin susceptibility, the existence of an enhanced density 
of states at the Fermi level while the Fermi surface sum 
rule includes the total number of tIc + 1 states. The work 
opens the way for further investigations as the systematic 
study of the effects of fluctuations in the pure Kondo 
lattice model or the incorporation of additional RKKY 
interactions in an extended Kondo lattice model. 
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